home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / dpbtrs.z / dpbtrs
Text File  |  1996-03-14  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. DDDDPPPPBBBBTTTTRRRRSSSS((((3333FFFF))))                                                          DDDDPPPPBBBBTTTTRRRRSSSS((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DPBTRS - solve a system of linear equations A*X = B with a symmetric
  10.      positive definite band matrix A using the Cholesky factorization A =
  11.      U**T*U or A = L*L**T computed by DPBTRF
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  14.      SUBROUTINE DPBTRS( UPLO, N, KD, NRHS, AB, LDAB, B, LDB, INFO )
  15.  
  16.          CHARACTER      UPLO
  17.  
  18.          INTEGER        INFO, KD, LDAB, LDB, N, NRHS
  19.  
  20.          DOUBLE         PRECISION AB( LDAB, * ), B( LDB, * )
  21.  
  22. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  23.      DPBTRS solves a system of linear equations A*X = B with a symmetric
  24.      positive definite band matrix A using the Cholesky factorization A =
  25.      U**T*U or A = L*L**T computed by DPBTRF.
  26.  
  27.  
  28. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  29.      UPLO    (input) CHARACTER*1
  30.              = 'U':  Upper triangular factor stored in AB;
  31.              = 'L':  Lower triangular factor stored in AB.
  32.  
  33.      N       (input) INTEGER
  34.              The order of the matrix A.  N >= 0.
  35.  
  36.      KD      (input) INTEGER
  37.              The number of superdiagonals of the matrix A if UPLO = 'U', or
  38.              the number of subdiagonals if UPLO = 'L'.  KD >= 0.
  39.  
  40.      NRHS    (input) INTEGER
  41.              The number of right hand sides, i.e., the number of columns of
  42.              the matrix B.  NRHS >= 0.
  43.  
  44.      AB      (input) DOUBLE PRECISION array, dimension (LDAB,N)
  45.              The triangular factor U or L from the Cholesky factorization A =
  46.              U**T*U or A = L*L**T of the band matrix A, stored in the first
  47.              KD+1 rows of the array.  The j-th column of U or L is stored in
  48.              the j-th column of the array AB as follows:  if UPLO ='U',
  49.              AB(kd+1+i-j,j) = U(i,j) for max(1,j-kd)<=i<=j; if UPLO ='L',
  50.              AB(1+i-j,j)    = L(i,j) for j<=i<=min(n,j+kd).
  51.  
  52.      LDAB    (input) INTEGER
  53.              The leading dimension of the array AB.  LDAB >= KD+1.
  54.  
  55.      B       (input/output) DOUBLE PRECISION array, dimension (LDB,NRHS)
  56.              On entry, the right hand side matrix B.  On exit, the solution
  57.              matrix X.
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. DDDDPPPPBBBBTTTTRRRRSSSS((((3333FFFF))))                                                          DDDDPPPPBBBBTTTTRRRRSSSS((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      LDB     (input) INTEGER
  75.              The leading dimension of the array B.  LDB >= max(1,N).
  76.  
  77.      INFO    (output) INTEGER
  78.              = 0:  successful exit
  79.              < 0:  if INFO = -i, the i-th argument had an illegal value
  80.  
  81.  
  82.  
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.